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ABSTRACT 

We report on observations of the superbubbles (SBs) Nil and N51D in the Large Magellanic Cloud 
(LMC) with Suzaku and XMM-Newton. The interior of both SBs exhibits diffuse X-ray emission, which 
is well represented by thin thermal plasma models with a temperature of 0.2-0.3 keV. The presence of 
nonthermal emission, claimed in previous works, is much less evident in our careful investigation. The 
3cr upper limits of 2-10 keV flux are 3.6 x 10" 14 ergs cm" 2 s" 1 and 4.7 x 10 -14 ergs cm" 2 s" 1 for Nil 
and N51D, respectively. The previous claims of the detection of nonthermal emission are probably due 
to the inaccurate estimation of the non X-ray background. We conclude that no credible nonthermal 
emission has been detected from the SBs in the LMC, with the exception of 30 Dor C. 
Subject headings: ISM: bubbles — ISM: individual (Nil, N51D) — X-rays: ISM 
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1. INTRODUCTION 

Massive stars are often found clustered in OB asso- 
ciations. Their high-speed stellar winds and the blast 
wave from core-collapsed supernova (SN) explosions col- 
lectively sweep up the ambient interstellar medium (ISM) 
to generate large (10-100-pc-scale) shell-like structures 
known as superbubbles (SBs) (e.g., McCray & Snow 
1979; Mac Low & McCray 1988). A tenuous cavity cre- 
ated inside an SB shell allows the blast wave of a sub- 
sequent supernova remnant (SNR) to propagate rapidly 
without deceleration for a long time. Furthermore, wind- 
wind and/or shock-cloud interactions, expected inside an 
SB, can maintain significant turbulence and magnetic in- 
homogeneities. Thus, SBs are believed to be very efficient 
site for the acceleration of relativistic cosmic rays (e.g., 
Bykov & Fleishman 1992; Parizot et al. 2004). 

The Large Magellanic Cloud (LMC) is particularly 
suitable for the stu dy of SBs, because of its low fore- 
ground extinction ([Dickey fc Lockmanl I1990T) and its 
known distance of ^50 kpc with a small uncertainty (e.g., 
Persson et al 2004). Us i ng RO SAT, Chu & Mac Low 
(1990) and [5 unne et aLl (|2001|) systematically studied 
several SBs in the LMC. They found that the cavities in- 
side the SBs exhibit diffuse X-ray emissions from rarefied 
hot gas with a temperature on the order of 0.1 keV. They 
also showed that the X-ray surface brightnesses are sig- 
nificantly higher than those expected from the standard 
pressure-driven SB model (Weaver et al. 1977). This 
suggests the presence of interior SNRs colliding with the 
inner walls of the SB shell. 

The detection of nonthermal X-rays has been rep orted 
in three LMC SBs, 3 Dor C (iBamba et al.|[200l. Ni l 
(|Maddox et al.l [20091) . and N51D ([Cooper et al.l l200l . 
30 Dor C exhibits the nonthermal emission with a clear 
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shell-like structure. The spectra obtained with Chan- 
dra and XMM-Newton were represe nted by a power-law 
with a photon index of T — 2-3 (Bamba et al. 2004), 
typical for syn chrotron X-rays from no nthermal SNRs. 
Using Suzaku, lYamaguchi et al.l (|2009| ) discovered the 
spectral cutoff in the power-law emission. They con- 
cluded that the synchrotron X-rays originate from the 
rapidly-expanding shell of the SNR with an age of (4- 
20) xlO 3 yr. 

On the other hand, th e validities of the oth er two re- 
sults may be doubtful. IMaddox et al.l ([2009D have re- 
cently observed Nil with Suzaku and have reported the 
detection of nonthermal X-rays above 2 keV. However, 
the claimed flux is comparable to that of 30 Dor C despite 
the fact that no nonthermal component w as detected 
by th e previous XMM-Newton observation ([Naze et al.l 
I2004D . Moreover, a morphology of the emission has been 
unclear. The Suzaku data should be therefore carefully 
reanalyzed. 

The nonthermal emission in N51D is claimed to be dis- 
tributed entirely in the SB cavity wi th a centrally peake d 
morphology in contrast to 30 Dor C ([Cooper et al.ll2004l) . 
The photon index of T ~ 1.3 is much lower than those 
expected for synchrotro n X-rays. On the oth er hand, in 
the preceding work by B omans et al.l ([2003D . the same 
XMM-Newton spectrum was well represented by a ther- 
mal plasma model without any additional nonthermal 
components. Thus, the existence of nonthermal emission 
in N51D is still controversial. Motivated by this issue, 
we have obtained a long-time exposure of this SB with 
Suzaku. 

The aim of this paper is to investigate for the exis- 
tence of diffuse hard X-ray emission. Therefore, other 
properties of SBs and their stellar contents are not fo- 
cused on. For the analyses, the Suzaku data are mainly 
utilized, since these provide a low and stable non- X-ray 
background (NXB), particularly for diffuse sources (Mit- 
suda et al. 2007; Tawa et al. 2008). In addition, we also 
refer to XMM-Newton archival data to compensate the 
limited angular resolution of Suzaku. 

First, we briefly introduce the characteristics of Nil 
( CTTTj) and N51D f fll~2"j) . In Section [2J we outline our 
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TABLE 1 

Log of observations appearing in this paper. 



Target 


Telescope/ 


Observation 


Aim point 


Observation 


Effective 


name 


Instrument 


ID 


a(J2000.0) 


<5(J2000.0) 


start date 


exposure (ks) 


Nil 


Suzaku/XIS 


501091010 


04 h 56 m 48 s 


-66°24'38" 


2006 Nov 7 


30.5 




XMM-Newton/EPIC 


0109260201 


04 h 55 m 58 s 


-66°25'53" 


2001 Sep 17 


32.5 (MOS)/25.1 (pn) 


N51D 


Suzaku/XIS 


804009010 


05 h 25 m 37 s 


-67°29'57" 


2009 May 17 


85.5/87.1/87.3 a 




XMM-Newton/EPIC 


0071940101 


05 h 26 m 03 s 


-67°28'52" 


2001 Oct 31 


30.5 (MOS)/23.1 (pn) 


NEP 


Suzaku/XIS 


500026010 


18 h ll m 18 s 


+66°00'45" 


2006 Feb 10 


88.5 



Note. — "The exposure times of the XISO/1/3. The XIS2 was out of operation during this observation. 



observation and data reduction procedure. The results 
for Nil and N51D are presented in Sections [3] and 01 
respectively. We discuss our results in Section [5j and 
give a summary in Section |6l The errors quoted in the 
text and tables refer to the 90% confidence intervals for 
one parameter, and the error bars in the figures are for 
la confidence, unless otherwise stated. 

1.1. Characteristics of Nil 

Nil (DEM L 34) , located at the northwest edge of the 
LMC (Henize 1956), con tains four OB associat ions, LH9, 
LH10, LH13, and LH14 (|Lucke &: Hodgdfl970l ). The Ha 
emission exhibits a complex structure consisting of sev- 
eral shells and filaments (e.g., Mac Low et al. 1998). The 
largest Ha shell with a diameter of ^120 pc is situated 
around the central cluster LH9. The age of LH9 is esti- 
mated to be ~3.5 Myr (Walborn et al. 1999), the oldest 
among the four clusters. It is suggested that the stellar 
activity of LH9 triggered a starburst at the periphery, 
leading to the formation of the other OB associations 
(e.g., Walborn & Parker 1992; Hatano et al. 2006). 

1.2. Characteristics of N51D 

N51D (DEM L 192) is located at the southeast corner 
of the supergian t shell LMC4. Two OB associations, 
LH51 and LH54 (jLucke fc Hodgdll970f ). are surrounded 
by a filamentary SB shell with a size of about 140 pc x 
120 pc. Both the OB clusters are estimated to have an 
age of -3 Myr (Oey & Smedley 1998). 

2. OBSERVATION AND DATA REDUCTION 

2.1. Suzaku 

The SBs Nil and N51D were observed with the Suzaku 
X-ray Imaging Spectrometer (XIS: Koyama et al. 2007) 
on 2006 November 7 and 2009 May 17, respectively. The 
observation IDs and aim points are summarized in Ta- 
ble [U The XIS consists of four X-ray charge coupled 
devices (CCDs). Three of them (XIS0, XIS2, and XIS3) 
are front-illuminated (FI) and the other (XIS1) is back- 
illuminated (BI). XIS0, 2, and 3 have a lower and more 
stable background level for extended emission, while the 
latter has superior sensitivity in the soft (<1.0 keV) 
X-ray band with a significantly improved energy reso- 
lution compared with previous X-ray imaging sensors. 
Combined with X-Ray Telescopes (XRTs; Serlemitsos et 
al. 2007), the field of view (FOV) of the XIS covers a 
~18'xl8' region with a half-power diameter (HPD) of 
~2'. During both the observations, the XIS was operated 
in the normal full-frame clocking mode with spaced-row 
charge injection technique. However, the XIS2 was out 
of operation during the observation of N51D due to dam- 
age, possibly caused by the impact of a micro meteorite. 



We reduced the data using the standard tools of 
HEADAS version 6.7. We first employed the revision 2.0 
and 2.4 data products for Nil and N51D, respectively. 
Both were reprocessed using the xispi software and the 
makepi files released on 2009 August 13. These files in- 
clude the latest calibration results for the charge transfer 
inefficiency and gain. We cleaned the reprocessed data 
in accordance with the standard screening criteri£0 and 
removed flickering pixels. We also excluded time inter- 
vals that suffered from telemetry saturation. After the 
filtering, the effective exposure times were obtained to 
be ~31 ks and ^86 ks for the observations of Nil and 
N51D, respectively. Only grade 0, 2, 3, 4, and 6 events 
were used in the following analysis. 

2.2. XMM-Newton 

Detailed information on the XMM-Newton observa- 
tions of both SBs is given in Table Q] The European 
Photon Imaging Camera (EPIC), aboard XMM-Newton, 
has two Metal Oxide Semi conductor (MOS) CCD arrays 
(Turner et al. 2001) and one pn CCD array (Struder et 
al. 2001). The latter has a higher quantum efficiency 
than the other CCD arrays. All three cameras have a 
larger FOV (^30'-diameter circle) than the XIS with 
a much better HPD of ~6" at the optical axis. This 
makes observations with the EPIC highly sensitive to 
point-like sources. During both the observations, the 
MOS and pn CCD arrays were operated in prime full 
window mode and extended prime full window mode, 
respectively. Thick filters were used in the observation 
of Nil, while thin ones were used in the observation of 
N51D. 

All the data were processed using version 9.0.0 of the 
XMM Science Analysis Software (SAS) with the lat- 
est calibration files. For the MOS, we selected X-ray 
events with patterns 0-12, which were passed through 
the #XMMEA_EM filter. For the pn, only events with 
patterns 0-4 and flag=0 were extracted. Light curves 
in the 10-15 keV band, where the effective area of the 
telescopes rapidly decreases, were accumulated from the 
whole FOV. To minimize any uncertainty due to the 
NXB, we removed the periods when the count rate was 
larger than 0.3 counts s _1 (MOS) and 1.0 counts s _1 
(pn) . The resultant exposure times are given in Table [TJ 

3. RESULTS FOR Nil 

3.1. Images and region selection 

Figures Ufa) and[llb) show XIS images of the Nil re- 
gion in the soft (0.5-2.0 keV) and hard (2.0-5.0 keV) 

5 http: / /heasarc. nasa. gov / docs /suzaku /processing/ criteria_xis . html 
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Fig. 1. — Images of Nil SB region, displayed on the J2000.0 coordinates. Panels (a) and (b) are Suzaku XIS images in the 0.5-2.0 keV 
and 2.0-5.0 keV bands, respectively. The latter is overlaid with contours of the former energy band smoothed with a Gaussian kernel of 
ct = 42". The ellipse and rectangle indicate the regions of source (SRC) and background (BGD) spectra, respectively. The small region 
shown as the gray circle in the panel (b) is excluded from the BGD spectrum. A panel (c) is XMM-Newton EPIC-MOS image in the 
2.0-5.0 keV band. The FOV of the Suzaku XIS and the regions for the spectral analysi s are shown. De t ected point sources are expressed 
as the small circles. Note that Shells 1, 2, and 3, labeled in the panel (a), are named by Mac Low ct al. (1998). 



X-ray bands, respectively. Data from the four detec- 
tors were combined. The images were binned by a factor 
of 8 and smoothed with a Gaussian kernel of a = 25". 
Extended emission can be clearly seen in the soft-band 
image, while no significant diffuse structure was found 
above 2.0 keV. 

As already mentioned in previous works (e.g., Naze et 
al. 2004), the soft emission peaks near the dense stellar 
cluster HD 32228 and is widely associated with the OB 
association LH9. This part is confined by the la rgest 
Ha shell of the Nil complex (|Mac Low et al.lll998l ). and 
thus it is considered that the soft X-rays originate from 
hot shocked gas within the SB shell that is blown by the 
cluster of stars. The emission is further extended to the 
north of LH9, the HII region N11B, and shells 2 and 3 
(see also Figure 3 of Mac Low et al. 1998). The X-rays 
from N11B are probably powered by stellar winds from 
the central young OB stars of LH10. Shells 2 and 3 lie 
to the north and northwest of N11B, respectively. These 
X-ray emissions with relatively low surface brightnesses 
a re confined by the fai nt Ha filaments. 

IMaddox et al.l (|2009[ ) claimed that nonthermal X-rays 
had been detected from the region around LH9. There- 
fore, we concentrate on the spectral analysis of the same 
region, an ellipse containing the entire SB around LH9 
with major and minor radii of 4.'5 and 3.'5, respectively — 
hereafter, "SRC" . The background (BGD) for the SRC 
spectra was taken from an off-source region, the rectangle 
shown in Figure [TJ but the CCD corners were excluded 
to remove X-rays from the calibration sources. 

A high-resolution EPIC image in the hard band is 
shown in Figure QJc). This image was made by merging 
two MOS images binned by a factor of 32 and smoothed 
with a Gaussian profile of a = 9."6. Similarly to the XIS 
image for the same energies (Figure [Hd), no extended 
emission was observed in the EPIC FOV. Point source ex- 
traction was performed using the ewavelet task of SAS 
with a detection threshold of 5a. The positions of the 
detected sources are indicated with small circles in the 
figure. Three discrete sources are located in the SRC re- 
gion, while one is located in the BGD region. The latter, 




0.5 1 2 5 

Energy (keV) 



Fig. 2. — Integrated EPIC spectra of three point sources included 
within the Nil SRC region. MOS1 (black), MOS2 (gray), and 
pn (black with circles) spectra are simultaneously fitted with an 
absorbed power-law model shown with solid and dashed lines. 

at (R.A., Decl.) = (4 h 57 m 28 s , -66°22'40"), is known to 
be associated with the cluster HD 268743, or Sk-66°41 
(N azc etaD[2004). Since this source is also clearly visible 
in the XIS image, we excluded it from the BGD region 
for spectral analysis. 

3.2. XMM-Newton spectra of point-like sources 

Here we analyze the EPIC data of the resolved point 
sources to take into account their contribution to the 
hard X-ray flux from the SRC region. Figure [2] shows the 
background-subtracted EPIC spectra, where the three 
point sources in the SRC region are merged to improve 
the statistics. The spectrum of each point source was 
extracted from the circular region with a radius of 15", 
while the background was taken from an annulus with 
inner and outer radii of 15" and 60", respectively. 

The spectra were phenomenologically fitted with an 
absorbed power-law model. The absorption cross sec- 
tions were taken from Morrison & McCammon (1983). 
The ele mental abundances w ere as sumed to be the val- 
ues m lAnders fc Grevessel (fl989l ). The fit was ac- 
ceptable (rVdof = 15/19) with N H = 2.9 (1.1-7.8) 
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TABLE 2 

2-10 KEV SURFACES BRIGHTNESSES OF Nil AND NEP REGIONS. 



Solid angle Flux [2-10 keV] Surface brightness [2-10 keV] 

(arcmin 2 ) (xlO - 13 ergs cm~ 2 s^ 1 ) (xl0~ 15 ergs cm~ 2 s" 1 arcmin" 2 ) 



SRC 


49.4 


2.2 (2.0-2.4) 


4.4 (4.0-4.8) 


(point sources) 




0.36 (0.31-0.41) 




(subtracted) 




1.8 (1.6-2.0) 


3.6 (3.2-4.1) 


BGD 


70.4 


2.8 (2.5-3.1) 


4.0 (3.6-4.4) 


NEP 


113.1 


5.0 (4.8-5.1) 


4.4 (4.2-4.5) 



Note. — The NXB is subtracted from each region, but the CXB is still included. 




0.5 1 2 5 10 0.5 1 2 5 10 

Energy (keV) Energy (keV) 



Fig. 3. — (a) XIS0 raw spectrum of the Nil SRC region (black) and NXB spectrum generated by xisnxbgen software (gray), (b) NXB- 
subtractcd photon spectra of the SRC (open circle) and BCD (filled square) regions, normalized by each solid angle. For comparison, the 
same spectrum from the NEP are shown with gray diamonds. The NEP spectrum was taken from an on-axis circular region with a radius 
of 6 arcmin. The dashed line indicates a power-law with photon index of 1.4, a typical value for the cosmic X-ray background. 



xlO 21 cm" 2 and T = 1.5 (1.1-2.2). The observed flux 
in the 2-10 keV band was obtained to be 3.6 (3.1-4.1) 
xl0~ 14 ergs cm~ 2 s^ 1 . 

3.3. Suzaku spectra of extended emission 

Since the emission from the SRC region is faint, back- 
ground reduction should be carefully performed to obtain 
an accurate result. We first subtract the NXB from the 
SRC and BGD raw spectra, and examine the difference 
in the hard X-ray fluxes between them (Section I3.3.1[) . 
Then, we subtract the other background components, 
such as the cosmic X-ray background (CXB) and local 
emission, and perform model fittings to the SRC spec- 
trum (Section [3.3.2[) . 

3.3.1. Search for hard X-rays 

Figure [3ja) shows the XIS0 raw spectrum of the SRC 
region. Instrumental emission lines of Mn-Ka, Ni-Ka 
and K/3, and Au-La can be observed at ~5.9 keV, 
~7.5 keV, ~8.3 keV, and ~9.7 keV, respectively. Using 
xisnxbgen software, we constructed the NXB spectrum 
from the night-Earth database. This was extracted from 
the same detector coordinates to the SRC region and 
sorted with the geomagnetic cut off rigidity for repro- 
ducibility. The intensities of the Ni-Ka lines were found 
to be 1.6 (1.1-2.1) xl(T 3 counts s" 1 and 1.6 (1.5-1.7) 
xl0~ 3 counts s _1 for the raw SRC and NXB spectra, 
respectively. Since these are highly consistent with each 
other, we can consider that the NXB spectrum was cor- 
rectly reproduced. The NXB for the BGD spectrum was 
also successfully constructed by the same method. 

After the NXB subtraction, the surface brightness 
spectra of the SRC and BGD regions were obtained as 




0.2 0.5 1 2 5 10 



Energy (keV) 

Fig. 4. — Effective-area ratios between the SRC and BGD regions, 
used for the vignetting correction. Solid black, solid gray, dashed 
black, and dashed gray represent XRT0, XRT1, XRT2, and XRT3, 
respectively. 

shown in Figure [3Jb) (shown as circles and squares, re- 
spectively). For comparison, we also present the spec- 
trum of another source-free region (shown as gray dia- 
monds) , taken from the observation of the north ecliptic 
pole (NEP). All the spectra are similar to each other at 
energies above ~3 keV, suggesting that the hard X-rays 
from the SRC region are dominated by the CXB. Tabled 
gives the 2-10 keV surface brightness of each region, de- 
termined using the data of all four XIS. No excess SRC 
flux was found. 

3.3.2. Spectral fitting 

The spectrum of the BGD region was subtracted from 
the SRC spectrum. However, since the vignetting of the 
Suzaku XRT significantly decreases the effective area at 
large off-axis positions (Serlemitsos et al. 2007), the BGD 
spectrum Ib(E) was corrected before the subtraction as 

I^E) =/ B (J5)~ •/(£), (1) 
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Fig. 5. — Background-subtracted XIS spectra of Nil, fitted with 
an APEC model (solid lines). Black and gray represent FI and BI, 
respectively. Contribution of the point sources determined from 
the EPIC spectra is given with the dotted lines. 



TABLE 3 

Best-fit spectral parameters for Nil. 



Parameter 

A^tj (cm" 2 ) a 

N™ c (cm- 2 ) 3 

kT e (keV) 

O (solar) 

Ne (solar) 

Mg (solar) 

Fe (solar) 

n e n p V (cm -3 ) 

-F£ bs (ergs cm -2 s -1 ) b 

L x (ergs s' 1 ^ 



Value 
4.3 X 10 2U (fixed) 

< 2.3 X 10 21 
0.30 (0.28-0.32) 
0.26 (0.15-0.47) 
0.39 (0.20-0.47) 
0.43 (0.22-0.89) 
0.12 (0.08-0.22) 
5.4 (2.6-19) x 10 58 
4.9 x 10~ 13 
2.2 x 10 35 



X 2 /dof 



99/105 = 0.94 



Note. "Absorption column in the Galaxy and LMC, 

where the solar and LMC abundances are respectively assumed. 
'Observed flux and luminosity in the 0.5-2.0 keV band. 



where fig and are the solid angles of the SRC 
and BGD regions. The correction factor f(E) is the 
energy-dependent effective-area ratio between the SRC 
and BGD regions. Using xissim software, we estimated 
the f(E) values for each XRT as given in Figure 21 This 
estimation also takes into account the contaminating ma- 
terial built up on the optical blocking filters of the XIS 
(Koyama et al. 2007). The differences in the results 
among the telescopes are mainly due to the scattering 
of the optical axes of the XRTs (see Figure 8 of Serlemit- 
sos et al. 2007). 

The BCD-subtracted SRC spectra are shown in Fig- 
ure The data from three FIs were combined be- 
cause their response characteristics are almost identi- 
cal. No significant signal was detected in the energy 
band above —3 keV, which is consistent with the re- 
sult in Section 13.3.11 K-shell emission lines of O VIII 
(-0.65 keV), Ne IX (-0.91 keV), and Mg XI (-1.34 keV) 
were clearly separate, indicating a thermal origin of the 
soft X-rays. Thus, we fitted the spectra with an optically 
thin thermal plasma model (APEC: Smith et al. 2001). 
The electron temperature (kT e ) and emission measure 
(n e n p V, where n e , n p , and V are the electron and pro- 
ton densities, and the emitting volume, respectively) 
were treated as free paramet ers. The elemental abun- 
dances relative to solar values (|Anders fe Grevesselll989|) 
were fixed to the mean LMC values of iRussell fe D opita 



(1992), but the Mg and Si abundances were assumed to 
be those of Hughes et a l. (1998) since the values from 
IRussell fe Dopital ((1992T) were more uncertain. The in- 
terstellar extinction in the Galaxy and LMC were sepa- 
rately considered. The Galactic absorption column den- 
sity with the solar abundances was fixed to be N§ — 
4.3 x 10 20 cm" 2 (jDickev fe Lockman|[T990T) . The other 
component (N^ MC ) was a free parameter, with the as- 
sumption of the LMC metal abundances. In addition, 
the point source model derived in Section [3.21 was given 
as a fixed component (with an independent absorption 
column). This model gave a best fit with kT e = 0.18 
(0.17-0.19) keV and x 2 /dof = 144/109. 

We next allowed the abundances of O, Ne, Mg, and 
Fe to vary freely (by using a VAPEC model). Then, 
we obtained a significantly improved fit with x 2 /dof 
= 99/105. The best-fit parameters are given in Ta- 
ble [3] The Fe abundance was found to be slightly lower 
than the mean LMC value (—0.3 solar: Russell & Do- 
pita 1992; Hughes et al. 1998). The observed flux of 
4.9 x 10~ 13 ergs cm" 2 s" 1 (in the 0.5-2.0 keV b and) 
is comparable with the r esult oflNaze eFal] (l200l but 
inconsistent with that of iMaddox et all (|2009f ) f— 2.5 x 



10~ i;i ergs cm~ 2 s _1 ), although the same data were an- 
alyzed. 

The derived electron temperature was slightly higher 
than the values in previous reports (—0.2 keV: Dunne 
et al. 2001; Naze et al. 2004; Maddox et al. 2009). We 
found, however, that the fit with a fixed kT e of 0.2 keV 
also yielded an acceptable x 2 /dof value of 111/106. This 
temperature difference does not severely affect the flux 
in the hard (>2 keV) X-ray band. 



4. RESULTS FOR N51D 



4.1. 



Images and region selection 

X-ray images of the N51D region are shown in Figure^ 
for the XIS soft X-ray band (a) , the XIS hard X-ray band 
(b), and the EPIC hard X-ray band (c). The binning 
factors and smoothing Gaussian kernels are same as those 
used in Figure [1] Similarly to Nil, the diffuse structure 
was only observed in the soft-band image. The bright soft 
emission is coincident with the OB association LH54 and 
is surrounded by the N51D Ha shell (see also Figure 1 of 
Cooper et al. 2004). 

The SRC regi o n was selected to be the same as that 
in ICooper et al.l ((2004) : an ellipse with major and mi- 
nor radii of 5.'8 and 4', respectively, and a 4'x4'-square, 
indicated in Figure [5] Although a small circular region 
around the Wolf-Rayet star HD 36402 was excluded from 
the spectral analysis in the previous work, here we in- 
cluded this region because no local structures were re- 
solved in the low-resolution XIS image. The BGD region 
was selected to be an off-source polygonal region, but the 
CCD corner irradiated by the calibration sources was ex- 
cluded. 

Using the EPIC image (Figure Wp), we searched for 
hard discrete sources. As a result, ten sources were de- 
tected in the SRC region, while no point source was found 
in the BGD region. These sources were not clearly re- 
solved in the XIS image. Therefore, we first analyzed the 
EPIC spectra of the point sources and then performed a 
model fit to the XIS data in the SRC region by adding 
the point-source component as a fixed model. 
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Fig. 7. — Same as Figure [2] but for N51D. Ten point sources 
within the SRC region are integrated. 



TABLE 4 

Flux distribution of the point sources within the N51D 
SRC region. 



Flux 


Number 


0.5-1.0 


2 


1.0-2.0 


3 


2.0-4.0 


4 


4.0-6.0 


1 



Note. — Observed fluxes are given in the spectral range of 
2-10 keV, and in unit of 10 -14 ergs cm -2 s _1 . 



4.2. XMM-Newton spectra of point-like sources 

EPIC spectra of the ten point sources in the SRC re- 
gion were extracted and integrated. The background 
spectra were extracted from an annular region around 
each point source. The resultant background-subtracted 
spectra are shown in Figure [71 The spectra were well 
reproduced with an absorption power-law model (x 2 /dof 
= 81/94). The best- fit absorption column and photon 
index were iV H = 1.3 (0.87-1.8) xlO 21 cur 2 and V = 1.7 
(1.6-1.8), respectively. The observed flux was obtained 
to be 2.4 (2.3-2.5) x 10" 13 ergs cur 2 s" 1 in the 2-10 keV 
band. The flux distribution of the ten sources is given in 
Table |H 



4.3. Suzaku spectra of extended emission 

The spectral analysis of the Suzaku data was performed 
by esse ntial ly the same procedure as that described in 
Section 13.31 An example of a raw spectrum in the SRC 
region and the corresponding NXB are compared in Fig- 
ure [Ha). Given that both spectra exhibit the same in- 
tensity of each fluorescent line, we considered that the 
NXB was successfully constructed. The NXB-subtracted 
spectra and the 2-10 keV surface brightnesses of the SRC 
and BCD regions were obtained as shown in Figure E^b) 
and Table El respectively. The SRC spectrum appears to 
show excess surface brightness in the hard X-ray band. 
However, upon subtracting the point-source flux from the 
SRC, we found that the remaining surface brightness is 
consistent with that in the BCD region. 

We applied the vignetting correction to the BCD spec- 
tra and then subtracted them from the SRC spectra. 
The correction factors f{E) were obtained as given in 
Figure |U The resultant SRC spectra are shown in Fig- 
ure [TO] The spectra of two active FIs were merged to 
improve the photon statistics. Several emission lines can 
be observed below ~1 keV, while the spectrum in the 
hard band is featureless. Therefore, we fitted the spec- 
trum with a VAPEC model plus the point-source compo- 
nent (a power-law) derived in Section 14.21 The elemen- 
tal abundances except for those of O, Ne, and Fe were 
fixed to the LMC values (Russell & Dopita 1992; Hughes 
et al. 1998). We assumed the Galactic absorption col- 
umn to be TVff = 6 . x 1 20 cm~ 2 , in accordance with 
iDickev fc~L ockmanl (|1990f) . The fit was acceptable with 
X 2 /dof = 92/103. The best-fit parameters and models 
are respectively shown in Table [6] and in Figure [10] with 
solid and dotted lines. Note that the point-source com- 
ponent was not fitted here. The model successfully repro- 
duced the spectral shape as well as the flux of the hard 
X-rays, and hence no additional nonthermal component 
w as required to impro ve the fit. 

ICooper et al.1 (]2004f ) argued that the EPIC spectrum 
of the diffuse emission exhibited thermal and nonthermal 
components, and the former accounted for ^70% of the 
total observed flux (in the 0.3-3.0 keV band) of 1.1 x 
10 -12 ergs cm -2 s _1 . With the given temperature and 
absorption column, the flux of the thermal emission was 
calculated to be ~6 x 10 -13 ergs cm~ 2 s^in the 0.5- 
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TABLE 5 

2-10 KEV SURFACES BRIGHTNESSES OF N51D REGION. 



Solid angle Flux [2-10 kcVJ 

(arcmin 2 ) (xl0 — 13 ergs cm -2 



Surface brightness [2-10 kcV| 
(xl0 — 15 ergs cm -2 s — 1 arcmin -2 ) 



5.2 (5.0-5.4) 
2.4 (2.3-2.5) 
2.8 (2.5-3.0) 
1.7 (1.5-1.9) 



- 1 ) 



5.9 (5.6-6.1) 

3.1 (2.9-3.4) 
3.5 (3.1-3.9) 



SRC 
(point sources) 
(subtracted) 

BCD 



48.4 




Fig. 9. — Same as Figure|9] but for the observation of N51D. 



TABLE 6 

Best-fit spectral parameters for N51D. 



Parameter 



value 



NX (cm" 2 ) 
AT™C (crn" 
kT e (keV) 
O (solar) 
Ne (solar) 
Fe (solar) 
n e n p V (cm -3 ) 
F£ bs (ergs cm -2 
Lx (ergs s' 1 ) 
X 2 /dof 



- 1 ) 



6.0 x 10 2U (fixed) 

< 5.2 x 10 20 
0.22 (0.20-0.23) 
0.14 (0.11-0.21) 
0.35 (0.26-0.49) 
0.18 (0.12-0.27) 
9.8 (7.0-14) x 10 58 
5.6 x 10" 13 
2.4 x 10 35 



92/103 = 0.89 



_ 


1 1 1 _ 

- | 


I 





0.5 



10 



Energy (keV) 



Fig. 10.— XIS spectra of N51D, fitted with an APEC model 
(solid lines). Black and gray represent FI and BI, respectively. 
The dotted lines represent the spectra of the point sources. 



2.0 keV band. This is comparable to our result. 

The enhancement of O a nd/or Ne abu n dance s was 
reported in t he w orks of iBomans et al.l (|2003l) and 
ICooper et al.l (|2004D . On the other hand, our spectrum 
does not show any evidence of a significant deviation 



from the typical LMC abundances, suggesting that the 
hot gas mainly originates from ISM with only slight en- 
richment by SN ejecta. 

5. DISCUSSION 

The X-ray spectra of the SBs Nil and N51D are well 
reproduced by single-component thin-thermal plasmas 
(with the models for the group of point sources). Fig- 
ure [11] shows the unfolded spectra of both SBs. The 
presence of additional nonthermal emission, argued in 
previous works (Nil: Maddox et al. 2009; N51D: Cooper 
et al. 2004), is much less evident i n our detail e d ana lyses, 
supporting the earlie r results by INaze et all (|2004l ) and 
IBomans et all (pOOl ). The 2-10 keV flux of the diffuse 
component Fd (± cd) for each SB can be calculated as 
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Fig. 11. — XIS-FI unfolded spectra (black crosses) of Nil (a) and N51D (b). The EPIC-pn spectra of the point sources are also shown 
with gray diamonds. Solid and da shed gray lines are th e mo del componen t s of t he diffuse and point sources, respectively. The spectra of 
"nonthermal emission" claimed by IMaddox et al. (2009) and Cooper ct al. (2004) (for Nil and N51D, respectively) are indicated with the 
black thick lines. The data points in the energies above ~ 2 keV are obviously far below the claimed fluxes although the point sources are 
not subtracted from the XIS spectra. 



where Pg, Fb, and Fps (and as, ob, ops) are the 2- 
10 keV fluxes (and uncertainties) of the SRC and BGD 
regions and that of the point sources. Using the values 
given in Tables [5] and [5j we obtain 3c upper limits of 
the fluxes to be 3.6 x 10 -14 ergs cm~ 2 s _1 and 4.7 x 
10~ 14 ergs cm" 2 s" 1 for Nil and N51D, respectively. 

For Nil, iMaddox et ail (|2009l) claimed that the hard 
X-rays, detected at energies of up to ~7 keV, required 
a power-law component with a photon index of T ~ 1.7 
and an observed flux of ^3.0 x 10~ 13 ergs cm -2 s _1 (in 
the 0.2-10 keV bandjS For N51D, the photon index and 
observed flux in the 0.3-3.0 keV band were reported to 
be r ~ 1.3 and ~3.3 x 10 -13 ergs cm -2 s _1 , respectively 
(|Cooper et al.ll200l . To directly compare these values 
with our results, these fluxes are modified to values in the 
2-10 keV band, giving 2.3 x 10~ 13 ergs cm" 2 s" 1 for Nil 
and 7.1 x 10~ 13 ergs cm~ 2 s _1 for N51D. These values 
are about one order higher than the upper limits that we 
have obtained. The spectral models for the claimed hard 
X-ray emissions are shown in Figure [TT] with black solid 
lines. Note that the point sources had been exclude d 
from the EPIC spectrum of N51D ([Cooper et al.|[200l . 
Nevertheless, both models give fluxes clearly exceeding 
those given by the XIS spectrum in which the point- 
source fluxes are contained. 

Given that the point-source flux in Nil is negligibly 
small (Fig. [Til) , the detection of hard X-rays claimed 
in the previous report was prob ably due to in c orrect 
background subtraction. In fact, IMaddox et al.l (2009) 
commented that the subtraction of Ni lines originating 
from the NXB was unsuccessful. The NXB spectrum of 
the XIS consists of emission-line and continuum compo- 
nents, and their in tensities are correlated with each other 
(jTawa et al.1 12008|) . Therefore, this also caused the in- 
complete subtraction of the continuum component in the 
previous reports. In our analyses, on the other hand, the 
fluorescent NXB lines are accurately subtracted (see Sec- 
tion 13.3. ip , and hence the background reduction is con- 
sidered to be successful. The problem concerning N51D 

6 The energy r a nge fr om which the flux was derived is not given 
in IMad dox ct al. (2009). Therefore, we independently estimated 
it from reported values of the absorption column density, photon 
index, and the ratio of the observed to unabsorbed fluxes. 



TABLE 7 

Comparison of nonthermal luminosities with other 
sources. 



Object 


L x (ergs s _1 ) a 


Reference 


Nil 


< 1 x 10 34 


This work. 


N51D 


< 1 x 10 34 


This work. 


30 Dor C 


2 x 10 35 


Yamaguchi et al. (2009) 


Arches 


1 x 10 34 


Tsujimoto et al. (2007) 


NGC 6334 


8 x 10 32 


Ezoe et al. (2006) 


RCW38 


5 x 10 32 


Wolk et al. (2002) 


Westerlund 1 


3 x 10 34 


Muno et al. (2006) 


Westerlund 2 


2 x 10 34 


Fujita et al. (2009) 


Note. — "Luminosity in the 2- 


-10 keV band. 



is analogous to the case of Nil, because the residual of 
the Al-Ka line due to i ncomplete backgrou nd subtrac- 
tion has been reported (jCooper et al.l I200"4l) . We point 
out that the analysis of faint extended sources, such as 
SBs in the LMC, is very sensitive to the estimation of 
the NXB and the vignetting effect, and thus should be 
performed with great care. 

Nonthermal X-ray emissions have so far been dis- 
covered from several Galactic star-forming regions, the 
Arches cluster (Tsujimoto et al. 2007), NGC 6334 (Ezoe 
et al. 2006), RCW 38 (Wolk et al. 2002), Westerlund 1 
(Muno et al. 2006), and Westerlund 2 (Fujita et al. 2009). 
In these cases, however, the emissions are correlated with 
relatively compact (1-10-pc-scale) regions around an OB 
association or in a molecular cloud core. 30 Dor C is 
another example of a star-forming region that exhibit s 
nonthermal (synchrotron) X-rays (|Bamba et al.ll200l . 
However, the X-ray shell can be interpreted as a single 
middle-aged SNR expan ding rapidly inside the SB cavity 
([Yamaguchi et al.l [20 09). Therefore, at present, there is 
no evidence for large-scale (~100 pc) nonthermal emis- 
sion distributed entirely in an SB. As given in Table [71 
the luminosities of the nonthermal X-ray emissions in 
Nil and N51D are, if any, less than those of 30 Dor C 
and some of the other star-forming regions. 

It is theoretically expected that efficient cosmic-ray 
acceleration can take place in SBs due to repeated 
SN shocks and/or magnetic turbulence (e.g., Bykov & 
Fleishman 1992; Parizot et al. 2004). Moreover, Butt 
& Bykov (2008) suggested that more than 10% of the 
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energy stored in an SB can be constantly transferred to 
accelerating cosmic-ray particles at the evolution stage 
of Nil and N51D (a few Myr). If this is the case, accel- 
erated electrons would quickly lose their energies via ra- 
diation so that no nonthermal X-rays would be detected. 
It is still an open issue whether or not SBs can efficiently 
accelerate cosmic-rays, but this is beyond the scope of 
this paper. Future detailed observations and the devel- 
opment of theoretical models are required to solve this 
problem. 

6. SUMMARY 

We have analyzed the Suzaku XIS data of the extended 
X-ray emission from the SBs Nil and N51D. The anal- 
ysis of the XMM-Newton EPIC data has also been per- 
formed to quantify the contribution from the point-like 
sources. In both SBs, the spectra of the diffuse com- 
ponents are well represented by thin thermal plasma at 
collisional ionization equilibrium. The electron temper- 
atures derived for Nil and N51D are ~0.3 keV and 
^0.2 keV, respectively. The elemental abundances are 
comparable with average LMC values, suggesting that 
the thermal emission mainly originates from the shock- 
heated ISM with only slight SN enrichment. 

Neither Nil nor N51D show any evidence for 
hard X-ray emission, contrary to the previous claims 



by iMaddox et ail (|2009D and iCooper et all (j2004l) . 
The 3cr upper limits of the 2-10 keV flux are 
3.6 x 10" 14 ergs cur 2 s" 1 for Nil and 4.7 x 
10 -14 ergs cm -2 s _1 for N51D. The published claims 
of the detection of nonthermal X-rays are thought to be 
simply due to inadequate background subtraction. Care- 
ful analyses (i.e., NXB and point-source subtraction, vi- 
gnetting correction, and so on) are required to obtain 
an accurate flux for faint extended sources, such as Nil 
and N51D. As the result of this work, no observational 
evidence for the nonthermal X-ray emission widely asso- 
ciated with SBs remains, with the exception of 30 Dor C 
(however, this is likely to be a single remnant of a re- 
cent SN explosion). It is still an open issue whether or 
not SBs can efficiently accelerate cosmic-rays. Observa- 
tions of more SBs with longer exposures as well as the 
development of theoretical calculations are required. 
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CXB. H.Y. is supported by the Special Postdoctoral Re- 
searchers Program in RIKEN, and the Grant-in-Aid for 
Young Scientists from the Ministry of Education, Cul- 
ture, Sports, Science and Technology (MEXT) of Japan. 
M.S. is a Research Fellow of Japan Society for the Pro- 
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